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Astrophysical Configurations with Background Cosmology: 
Probing Dark Energy at Astrophysical Scales 



o 
o 

(N 

o 
O 



6 



(N 
> 
O 
00 

a^ 

00 

O 
l> 
O 



A. Balaguera-Antolmez^*, D. F. Mota^f and M. Nowakowski^J 

^ Max Planck Institute fiir Extraterrestrische Physik, D-85748, Garching, Germany 

Institute for Theoretical Physics, University of Heidelberg, 69120 Heidelberg, Germany 
^ Departamento de Ftsica, Universidad de las Andes, A. A. 4976, Bogota, D.C., Colombia 



1 February 2008 



ABSTRACT 

We explore the effects of a positive cosmological constant on astrophysical and cos- 
mological configurations described by a polytropic equation of state. We derive the 
conditions for equilibrium and stability of such configurations and consider some astro- 
physical examples where our analysis may be relevant. We show that in the presence 
of the cosmological constant the isothermal sphere is not a viable astrophysical model 
since the density in this model does not go asymptotically to zero. The cosmological 
constant implies that, for polytropic index smaller than five, the central density has 
to exceed a certain minimal value in terms of the vacuum density in order to guar- 
antee the existence of a finite size object. We examine such configurations together 
with effects of A in other exotic possibilities, such as neutrino and boson stars, and 
we compare our results to N-body simulations. The astrophysical properties and con- 
figurations found in this article are specific features resulting from the existence of a 
dark energy component. Hence, if found in nature would be an independent probe of 
a cosmological constant, complementary to other observations. 
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1 INTRODUCTION 

Some of the most relevant propert ies of the un iverse have been established through astronomical data associated to light 
curves of distant Sup ernova la ( Rie ss et al. 2004 1 . the temperature anisotropics in the cosmic microwa ve background radiation 
i Spergel et al. 20061 ) and the matter power spectrum of large scale structures ( Tegmark et al. 2004h . Such observations give 
a strong evidence that the geometry of the universe is flat and that our Universe is undertaking an accelerated expansion at 
the present epoch. This acceleration is attributed to a dominant dark energy component, whose most popular candidate is 
the cosmological constant, A. 

The present days dominance of dark energy make us wonder if this component may affect the formation and sta- 
bility of large astrophysical structures, whose phy sics is basically Newtonian. This is in fact an old question put for 
ward already by Einstein ( Einstein fc Straus 1945h and p ursued by many other authors ( Noerdlinger fc Petrosian 1971 



fernin Nagirner & Starikova 2003'; 'Manera & Mota 2005'; ^Nu nes fc Mota 20061 : iBarvshev.Chernin fc Teerikorpi 20011 ) and 
agramanova.Kunz fc Lacmmcrzahl 2006; Jotzcr & Serena 200(^). In general, the problem is rooted in the question whether 
the expansion of the universe, which in the Newtonian sense could be understood as a repulsive force, affects local astrophys- 
ical properties of large objects. The answer is certainly affirmative if part of the terms responsible for the Universe expansion 
survives the Newtonian limit of the Einstein equations. This is indeed the case of A which is part of the Einstein tensor. In 
fact, as explained in the main text below, all the effects of the universe expansion can be taken into account, regardless of the 
model, by generalizing the Newtonian limit. This approach allow us to calculate the impact of a given cosmological model on 
astrophysical structures. 
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Althoush, there are s everal candidates to dark energy which have th eir own cosmoloRical signature, e.g. jKoivisto fc Mota 2007 ; 
Daly fc Diorgovski 2004IWang fc Mukheriee 20041: iKoivisto fc Mota'2 006: K oivisto fc Mota 2007al: iMota fc Shaw 20071 ) and 



i Seo fc Eisenstein 200^ : 1 Brookfield et al. 20061 : 



Daly fc Diorgovski 20031 : iKoivisto fc Mota 2007bl : lMota fc Shaw 2006[ rin this 



paper we will investigate the ACDM model only. Such consideration is in fact not restrictive and our results will be common 
to most dark energy models. At astrophysical scales and within the Newtonian limit one does not expect to find impor- 
tant differences among the different dark energy models. This, however should not be interpreted as if A has no effect 
at smaller astrophysical scales. In fact, the effects of a cosmological constant on the equilibrium and stability of astro- 
physical structures is not negligible, and can be of relevance to describe features of astrophysical systems such as globular 
clusters, galaxy clusters or even galaxies |Ch«rnii^_et^L_^007 ^ lorio 20051 : Balaguera-Antolmez.Bohmer fc Nowakowski 2006b ; 
Nowakowski fc Sanabria fc Garcia 20021 : ICardoso''fc''Guaitieri 2006l ) . Motivated by this, we investigate the effects of a dark 



energy component on the Newtonian limit of Einstein gravity and its consequences at astrophysical scales. 

In this article we investigate how the cosmological constant changes certain aspects of astrophysical hydrostatic equilib- 
rium. In particular we search for specific imprints which are unique to the existence of a dark energy fiuid. For instance, the 
instability of previously viable astrophysical models when A is included. We explore such possibility using spherical config- 
urations described by a polytropic equation of state (e.o.s) p ~ p''. The polytropic equation of state derives its importance 
from i ts success and consiste ncy, and it i s widely used in determini ng the properties of gravitational structures ranging from 



stars 



Chandrasekhar 19671 ) to galaxies (jBinnev fc Tremaine 19871 ). It leads to an acceptable description of the behavior of 



astrophysical objects in accordance with observations and numerical simulations jKennedv fc Bludman 199S Gruzinov 2000l : 



Kaniadakis.Lavagno fc Quarati 19961 : ISadeth fc Rephaeh 20041 : Ipinzon fc Calvo-Mozo 200ll:lR uffet et al. 19961 ). The descrip 



tion of such configurations can be verified in the general relativistic f ramework (Herrera fc Barreto 2 0031 ) and applications of 
these models to the dark energy problem have in fact been explored I Mukhopadhvav fc Ray 2005t ) . 

The effect of a positive cosmological constant can be best visualized as a repulsive non local force acting on the mat- 
ter distribution. It is clear that this extra force will result into a minimum density (either central or average) which is 
possible for the distribution to be in equilibrium. This minimum density i s a crucial crossin g point: below this value no 
matter can be in equilibrium, above this value low density objects exist jLahav et al. 199ll ). Both effects are novel fea- 
tures due to A. We will demonstrate s uch in equalities, which are generalizations of correspond ing inequalities found in 
( Nowakowski fc Sanabria fc Garcia 2002h and ( Balaeuera-Antoh'nez.Bohmer fc Nowakowski 2005al ). for every polytropic in- 
dex n. However, the most drastic effect can be found in the limiting case of the polytropic equation of state, i.e, the isothermal 
sphere where the polytropic index n goes to infinity. This case captures, as far as the effects of A are concerned, many 
feat ures also for higher, but finite n . The model of the isothermal sphere is often used to model galajcies and galactic clus- 
ter s ('Natar aian fc Lvnden-Bell 1997 ) and used in describing effects of gravitational lensing i Kawano et al. 2004 : Sereno 20051 : 
[Macc io 20oir Herein lies the importance of the model. Regarding the isothermal sphere we will show that A renders the 
model unacceptable on general grounds. This essentially means that the model does not even have an appealing asymptotic 
behavior for large radii and any attempt to definite a physically acceptable radius has its severe drawbacks. 

The positive cosmological constant offers, however, yet another unique opportunity, namely the possible existence of 
young low density virialized objects, understood as configurations that have reached virial equilibrium just at the vacuum 
dominated epoch (in contrast to the structures forming during the matter dominated era, where the criteria for virialization 
is roughly p ~ 200pcrit). This low density hydrostatic/virialized objects can be explained again due to A which now partly 
plays the role of the outward pressure. 

The applicability of fluid models, virial theorem and hydrostatic equilibrium to lar ge astrophysical bodies has been dis- 

cussed many times in the literature. For a small survey on this topic we refer the reader to (|jackson 197GI : Balaeuera-AntoKnez.Mota fc Nowako' 
where one can also find the relevant references. 

It is interesting to notice that dark matter halos represent a constant density background which, in the Newtonian limit, 
objects embedded in th em feel the analog to a negative c osmological constant. The equilibrium analysis for such configurations 
has been performed in ( Umemura fc Ikeuchi 198a : Horedt 20 00). A negative A will just enhance the attractive gravity effect, 
whereas a positive one opposes this attraction. As a result the case A > reveals different physical concepts as discussed in 
this paper. 

The article is organized as follows. In the next section we introduce the equations relevant for astrophysical systems as 
a result of the weak field limit and the non-relativistic limit of Einstein field equations taking into account a cosmological 
constant. There we derive such limit taking into account the background expansion independently of the dark energy model. 
In section 3 we derive the equations governing polytropic configurations, the equilibrium conditions and stability criteria. In 
section 4 we describe the isothermal sphere and investigate its applicability in the presence of A. In section 5 we explore some 
examples of astrophysical configurations where the cosmological constant may play a relevant role. In particular, we probe 
into low density objects, fermion (neutrino) stars and boson stars. Finally we perform an important comparison between 
polytropic configurations with A and parameterized density of Dark Matter Halos. We end with conclusions. We use units 
Gjv = c = 1 except in section 5.4 where we restore Gn and use natural units ?!, = c = 1. 
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2 LOCAL DYNAMICS IN THE COSMOLOGY BACKGROUND 



The dynamics of the isotropic and homogeneous cosmological background is determined by the evolution of the (dimensionless) 
scale factor given through the Friedman-Robertson- Walker line element as solution of Einstein field equations, 



a{t) 



--7T[p{t) + 3pit)], 



a{t) 



Hit)' 



(1) 



corresponding to the Raychaudhury equation and Friedman equation, respectively. The total energy density p is a contribution 
from a matter component - baryonic plus dark matter - (pmat ~ radiation (prad ~ and a dark energy component 

(px ~ a"-^'"' with p = ujxp)- The function /(a) is given as 



-da , 



(2) 



where the term tJx(a) represents the equation of state for the dark energy component. The case Wx = — 1 corresponds to the cos 
mological constant px — Pvac ~ A/Svr. The effects of the background on virialized structures ca n be explored through the New 



tonian limit of fiel d equations from which one can derive a modified Poisson's equation (see e.g. ([Noerdlinger fc Petrosian 1971 



iNowakowsk i 2001 1). Recalling that pressure is also a source for gravity, the gravitational potential produced by an overdensity 
is given by 



V^$ = 47r(pt + 3Pt), 



(3) 



where pt = 5p + p and Pt — SP + P. Where 5p is the local overdensity with respect to the background density p. Notice that 
equation Q reduces to the usual Poisson equation, V^"I> = 4Tv5p, when non relativistic matter dominates the Universe, and 
Sp ^ p. However, at present times, when dark energy dominates, the pressure is non-negligib l e and 5P might even be no n 
zero, such as in the case of quintessence models IjMaor fc Lahav 20051 : IWang fc Steinhardt 19981 : iMota fc van de Bruck 2004 ). 
In this work, however, we will focus in the case of an homogeneous dark energy component where 5p = SP = 0. With this in 
mind, one can then write the modified Poisson equation as 



4.Sp^3% 
a{t) 



(4) 



Note that this equation allows one to probe local effects of different Dark Energy models through the term a /a given in Eq.([T]). 

Since we will be investigating the configuration and stability of astrophysical objects nowadays, when dark energy domi- 
nates, it is more instructive to write the above equations in terms of an effective vacuum density i.e. 



STrPvacI 



where by using ((T} p! 

^cdm 



a), 



vac ' 



a) has been defined as 



off , 

Pv!^c(a 



V fivac / ' 



+ (1 + ScJx 



-/(a) 



(5) 



(6) 



^^vac 

which reduces to pvac for uj^ = —1 and negligible contribution from the cold dark matter component with respect to the over- 
density 5p. With $grav being the solution associated to the pure gravitational interaction, the full solution for the potential 
can be simply written as 



$(r,a) = $g 



/ \ 4 off / ^ 2 

An - gl-p^acla)'^ , 



5p{v') 
|r — r'l 



(7) 



which defines the Newton-Hooke space-time for a scale factor close to the present time (vacuum dominated epoch), uj^ = —1 
and ricdm < fivac ( Gibbons fc Patricot 20031 : 1 Aldrovandi et al. 19981 ). For a ACDM universe with ficdm = 0.27 and f^vac = 0.73 
we get p°f(,(today) = 0.81pvac: that is, the positive density of matter which has an attractive effect opposing the repulsive one 
of A reduces effectively the strength of the 'external force' in ([5]). Note that, although in the text we will use the notation 
Pvac which would be valid in the case of a Newton-Hook space time, it must be understood that we can replace pvac by 
Pvfc(today) = O.Slpvac for a ACDM background. 

Given the potential $(r, a), we can write the Euler's equation for a self-gravitating configuration as 

d{vi) 



dt 



+ d^p + pd,$ = 



, where (vi) is the (statistical) mean velocity and p is the total energy density in the system. We can go beyond Euler's equation 
and write down the tensor virial equation which red uces to its scalar version for spherical configurations. Th e (scalar) virial 
equation with the background contribution reads as (jBalaguera-Antolmez fc Nowakowski 20051 : ICaimmi 20071 ) 



= 2r + 



3n + |7rp°fc(a)J - 



p(r ■ 



dA, 



(8) 
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where W^'^^^ is the gravitational potential energy defined by 



p(r)*grav(r)d r, 



(9) 



is the moment 



and T = ^ p{v'^)d'^T is the contribution of ordered motions to the kinetic energy. Also, X = pr~a"r 
of inertia about the center of the configuration and 11 = pd'^r is the trace of the dispersion tensor. The full descrip- 
tion of a self gravitating configuration is completed with an equation for mass conservation, energy conser vation and an 
equation of state p = p(Sp,s). If we ass ume equilibrium via X ~ 0, we obtain the known virial theorem i Jackson 1970l : 
Balaeuera-Antolmez fc Nowakowski 20051 ) 



2r + an 



(10) 



where we have neglected the surface term in ((Sjl , which is valid in the case when we define the boundary of the configuration 
where p = 0. 

With given in one must be aware that an equilibrium configuration is at the most a dynamical one. This 
is to say that the 'external repulsive force' in ([§} is time dependent through the inclusion of the background expansion 
and so arc the term s in Th is leads to a viol ation of energy cons ervation, which also occu rs in the virialization process 
( Wang 2006: Maor fc Lahav 200 5: Caimmi 2007: Mota fc van de Bruc k 2004; Shaw fc Mota 2007A '). Traced over cosmological 



times, this implies that if we insist on the second derivative of the inertial tensor to be zero, then the internal properties of 
the object like angular velocity or the internal mean velocity of the components will change with time. Even if in the simplest 
case, one can assume that the objects shape and its density remain constant. Hence, equilibrium here can be thought of as 
represented by long time averag es in which case the second derivative al s o vanishes, not because of constant volume an d 
density, but because of stability ( Balaguera- Antolmez fc No wakowski 20061 : Balaguera-Antoh'nez.Mota fc Nowakowski 20061 ). 

The expressions derived in the last section, especially ([5} and (|10[l can be used for testing dark energy models on 
configurations in a dynamical state of equilibrium. However, in these cases, one should point out that, in this approach 
there is no energy conservation within the overdensity: dark energy flows in and out of the overdensity. Such feature is a 
consequence of the assumption that dark energ y does not cluster at small scales (hom ogeneity of dark energy). This is in fact 
the most common assumption in the literaturedWang fc Steinhardt 199S : Chen fc Ratra 2004 : Horellou fc Berge 20051 ), with 
a few exceptions investigated in (|Wang 200d : lMaor fc Lahav 2005 : I Caimmi 20071 : I Mot a fc van de Bruck 20ol v 

In this paper, we will concentrate on the possible effects of a background dominated by a dark energy component 
represented by the cosmological constant at late times {z ^1). It implies that the total density involved in the definitions of 
the integral quantities appearing in the virial equation can be approximated to p ~ Sp. In that case the Poisson equation 
reduces to the form V^'I' = AnSp — STrpvac. As mentioned before, the symbol pvac has no multiplicative factors in the case of 
a Newton-Hooke space time, while for a ACDM model it must be understood as pvac ^ O.Spvac- As the reader will see, the 
most relevant quantities derived here come in forms of ratio of a characterizing density and pvac, and hence the extra factor 
appearing in the ACDM can be re-int roduced in the characte rizing density. For general consideration of equilibrium in the 
spherical case see IIBohmer 2004) and (Bohmer fc Harko 20051 ). while the quasi spherical collapse with cosmological constant 
has been discussed in IjPebnath et al 20061 ). 



3 POLYTROPIC CONFIGURATIONS AND THE A-LANE-EMDEN EQUATION 

We can determine the relevant features of astrophysical systems by solving the dynamical equations describing a self gravitating 
configuration (Euler's equation, Poisson's equation, continuity equations). In order to achieve this goal we must first know the 
potential $ to be able to calculate the gravitational potential energy. To obtain $, one must supply the density profile and 
solve Poisson's equation. In certain cases, the potential is given and we therefore can solve for the density profile in a simple 
way. Here we face the situation where no information on the potential (aside from its boundary conditions) is a vailable and 



we also do not have an apriori information about the density profile (see for instance (jBinnev fc Tremaine 19871 ) for related 
examples). In order to determine both, the potential and the density profile, a complete description of astrophysical systems 
required. This means we need to know an equation of state p — p{p) (here we change notation and we call p the proper density 
of the system). The equation of state can take several forms and the most widely used one is the so-called polytropic equation 
of state, expressed as 

p=Kp^, 7EEl + i, (11) 

n 

where 7 is the polytropic index and k is a parameter that depends on the polytropic index, central density, the mass and the 
radius of the system. The exponent 7 is defined as 7 = (cp— c)/(cv— c) and is associated with processes with constant (non-zero) 
specific heat c. It reduces to the adiabatic exponent if c = 0. The polytropic equation of state was introduced to model fully 
convective configurations. From a statistical point of view, Eq. (|ll|l represents a coUisionless system whose distribution function 
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Figure 1. Solutions of ALE equation for different (c and different index n ranging from n = 3 (red, short-dashed hne), n = 4 (blue, long- 
dashed line) and n = 5 (green, dot short-dashed line). 



can be written in the form / = f{E) ~ ^"~3/2^ with E = (p— (l/2)mv^ being th e relative energy and (/)(r) = $o — being 
the relative potential (where $o is a constant chosen such that <j){r = R) = 0) (|Binnev fc Tremaine 1987^ . In astrophysical 
contexts, the polytropic equ ation of state is widely used to describe astrophysical systems such as the su n, compact objects, 
galaxies and galaxy clusters ( Chandrasekhar 1967 : Binnev fc Tremaine 1987 : Kennedv fc Bludman 19991 ) ■ We now derive the 
well known Lane-Emden equation. We start from Poisson equation and Euler equation for spherically symmetric configurations, 
written as 



d2$ 2 d$ 

-^TT ^ r' ^ 47rp — 87r/9vac, 

dr^ r ar 



dp d$ 
dr ^ dr 



(12) 



This set of equations together with Eq. (|lip can be integrated in order to solve for the density in terms of the potential as 

7-1 



p(r) 



Pc 



1 



K7 



1-7 



mr) - $(0)) 



(13) 



where pc is the central density. In view of eq. (|13p in conjunction with eq. ((7} it is clear that pvac will have the effect to 
increase the value of p{r). Therefore the boundary of the configuration will be located in a greater R as compared to the case 
Pvac = 0. In order to determine the behavior of the density profile, we again combine Eq (|12|l and Eq. (|lip in order to eliminate 
the potential <l>(r). We obtain 



n y p J K(n + 1) 

where we defined the function 



C = C(r) 



Pvac 

p{r) 



(14) 



(15) 
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We can rewrite Eq (|14|l by introducing the variable ip defined by p = Pc4>", where pc is the central density. We also introduce 
the variable ^ = r/a where 

(16) 




is a length scale. Ea. (|14|l is finally written as (|Balaguera-Antolmez.Bohmer fc Nowakowski 2005al : I Chandrasekhar 19671 ) 



This is the A-Lane-Emden equation (ALE) . Note that for constant density, we recover p = 2pvac as the first non trivial solution 
of ALE equatio n. This is consistent with the results from virial theorem for constant density spherical objects which tell us 



that p ^ 2pvac (jNowakowski fc Sanabria fc Garcia 20021 ). Note that using Eq. (|13p we can write the solution V'(C) with the 
explicit contribution of pvac as 

= r/a) = 1 - (47raVc)"' ($grav(r) - $g.av(0)) + 6Ccf, (18) 

so that for a given r smaller than the radius we will obtain 

^/'(r/a) > V(r/a)A=o. (19) 

as already pointed out before. Then the differential equation (|17p must be sol ved with the initial conditions tp(0) = 1, i/^ 'jO) = 



0, satisfied by (|18p . Numerical solutions were obtained for the first time in (jBalaguera-Antohnez fc Nowakowski 2005f ). The 
solutions presented in Fig. [1] are given in terms of the ratio p/pvac for n = 3,4 and n = 5. This choice of variables are useful 
also since pvac sets a fundamental scale of density (the choice po = pvac will be explored for the isothermal sphere, where figure 
[T]will be helpful for discussions). The radius of a polytropic configuration is determined as the value of r when the density of 
matter with the e.o.s pip vanishes. This happens at a radius located at 

R = a^i such that ip{^i) = 0. (20) 

Note that equation (|17p yields a transcendental equation to determine ^i. Also one notes from Fig. [T] that not all values of 
yield allowed configurations in the sense that we cannot find a value of such that i/'(5i) = 0. This might not be surprising 
since for n = 5 and A = we find the situation where the asymptotic behavior is p — > as ^ — + oo (we consider this still as an 
acceptable behavior). There is, however, one crucial difference when we switch on a non-zero A. For A 7^ not only we cannot 
reach a definite radius but the derivative of the density changes sign and hence becomes non-physical. The situation for the 
cases n 55 5 is somewhat similar to the extreme case of n — > 00 (isothermal sphere). Clearly, these features are responsible of 
the last term in Eq. (|18p . which for high values of ( may become dominant over the remaining (gravitational) terms. We will 
discuss this case in section four where we will attempt another definition of a finite radius with the constraint V'(C)' < 0- For 
now it is sufficient to mention that, as expected, the radius of the allowed configurations are larger than the corresponding 
radius when A = 0. 



3.1 Equilibrium and stability for polytropes 

In this section we will derive the equilibrium conditions for polytropic configurations in the presence of a positive cosmological 
constant. We will use the results of last section in order to write down the virial theorem. The total mass of the configuration 
can be determined as usual with M = J pd'^r together with Ea. (|17p . One then has a relation between the mass, the radius 
and the central density: 

R = M^/Vc-'/Vo(Cc; n) = (Mri)i/3(4^)i/\y Vo(Cc; n) (21) 
where 

. 1 , t ' 
51 I ( / ^2 



/o(Cc;n) = (^^j i^j erm^j , (22) 

Note that we have introduced the cosmological constant in the equation for the radius, leading to the appearance of the 
astrophysical length scale (Mr^)^''^ (with ta = A'^''^ = (87rpvac)"^/^ = 2.4 x 10^ (!r2vac/!.^)"^''^Mpc ^ 4.14 x 10^ Mpc for the 
concordance values 57vac =0.7 and h — 0.7). This scale has been already found in the context of Schwarzschild - de Sitter metric 
where it is the maximum allowed radius for bound orbits. At the same time it is the scale of the majdmum radius for a self grav - 
itating spherical and homogeneous configuration in the presence of a positive A ( Balaeuera- Antolmez fc Nowakowski 2005l l. 



This also let us relate the mean density of the configuration with its central density and/or cosmological parameters as 

P= (3/47r/(?)pc = (3/27rCc/(?)Pvac = (3!:!vac/27rCc/o )Pcrit. 
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Similarly we can determine the other relevant quantities appearing in the equations for the energy and the scalar virial 
theorem (|10p . For the traces of the moment of inertia tensor and the dispersion tensor 11 we can write 

J=Mi?Vi, n = Acpc"M/2, (23) 

where the functions /i,2 have been defined as 



/i(Cc;n) 



/2(Cc;n) 



■Ci(Cc 



(24) 



/i(Cc), such that for a given mass 



using l|2ip . These functions are numerically determined in the sequence ip{^\ Cc) 
we obtain the radius as i? = a(Af;^i)5i. 

Let us consider the virial theorem pO[) for a polyt ropic co nfigurati on. The gravitational potential energy W^rav 
obtained following the same arguments shown in (| Chandrasekhar 19671 ). The method consist in integrating Euler's equation 



and solve for $grav, then using Eq.(|9} one obtains W^"^"^. The final result is written as 



hp 



i(n + l)n + Ivrpvac (X - MR") , 



(25) 



2R 2' ' 3 

To show the behavior of W^"^^^ with respect to the index n, we can solve the the virial theorem (|10|l for 11 and replace it in 
Eq.((25l). We obtain 

(26) 



5 — n 



Pvac 



-(5 + 2n)/i 



This expression shows the typical behavior of a n = 5 polytrope (even if pvac 7^ 0): the configuration has an infinite potential 
energy, due to the fact that the matter is distributed in a infinite volume. The energy of the configuration in terms of the 
polytropic index can be easily obtained by using (|10p . (|25|l and (I26[l 



E = W^'^" + lyrpvacX + nn = - 



3 - n \ 
5-n) ~R" 



(n) 



(27) 



One is tempted to use i? < as the condition to be fulfilled for a gravitationally bounded system. For pvac = we recover 
the condition 7 > 4/3 (n < 3) for gravitationally bounded configurations in equilibrium. On the other hand, for pvac 7^ 
this condition might not be completely true due to the following reasoning: The two-body effective potential in the presence 
of a positive cosmological constant does not go asymptotically to zero for large distances, which is to say that _E < is not 
stringent enough to guarantee a bound system. Therefore, we rather rely on the numerical solutions from which, for every n, 
we infer the value of An such that 



Pc ^ An pv 



(28) 



This gives us the lowest possible central density in terms of pvac- The behavior of the ("ciit, the functions fi, the solution 
^i{(c = Ccrit, n) and the values of An are shown in figlH Note that this inequality can be understood as a generalization of the 
equilibriu m condition p > ^pvac, which, when applied for a spherical homogeneous configurations yields A = 2 with Q — p — 



constant ( Balaeuera-Antolfnez fc Nowakowski 20051 : iBalaguera-Antolmez.Bohmer fc Nowakowski 2005a 1. 

Note that, at n = 5 the radius of the configura tion becomes undefined, as well as the energy. A n — 5 polytrope is highly 



concentrated at the center 



Chandrasekhar 19671 1. No criteria can be written since even for pv 



there is not a finite 



radius. But it is this high concentration at the center and a smooth asymptotic behavior which makes this case still a viable 
phenomenological model if A is zero. On the contrary for non-zero, positive A the solutions start oscillating around 2pvac 
which makes the definition of the radius more problematic. For n —> 00 the polytropic e.o.s describes an ideal gas (isothermal 
sphere) . Since in this limit the expressions derived before are not well definite, this case will be explored in more detail in the 
next section. In spite of the mathematical differences, the isothermal sphere bears many similarities to the cases n ^ 5 and 
our conclusions regarding the definition of a radius in the n —> 00 case equally apply to finite n bigger than 5. 



3.2 Effects with generalized dark energy equation of state 

In the last section we have explored the effects of a dark energy-dominated background with the equation of state ljx ~ — 1. 
Other dark energy models are often used with uj^ = —1/3 and lj^ = —2/3 or even uj^ < —1, in the so-called phantom regime, 
or even a time dependent dark energy model (quintessence). A simple generalization to such models can be easily done by 
making the following replacement in our equations 

C^aaU = -^Cv{a)a-f'-^\ (29) 
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Figure 3. Density of a polytropic configuration for ijc = 10 ^ and different polytropic index in a modified Newton-Hook space time 
with a dark energy equation of state oJx = —2/3 (black, dots) and lo^ = —2 (red, short dash). Compare with Fig[T] 
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where 77(a) = 1 + 3u>x{a), and where the function /(a) is defined in Eq.©. Note that for this generalization to be coherent 
with the derivation of ALE in (|17|) . one needs to consider that the equation of state is close to —1, so that there is almost 
no time-dependence, and the energy density for dark energy is almost constant. This is indeed the case for most popular 
candidates of dark energy specially at low redshifts 2 < 1. Also, notice once again, that we are still assuming an homogeneous 
dark energy component which flows freely to and from the overdensity. Hence, violating energy conservation inside it. Clearly, 
other models of dark energy will posses dynamical properties that the cosmological cons tant does not have. For instance, we 
could allow some fraction o f dark energy to take part in the collapse and virialization ( Maor fc Lahav 20051 : Caimmi 2007 : 



Mota fc van de Bruck 2004 '). which would lead to the presence of self and cross interaction terms for dark energy and the 



(polytropic like) matter in Euler equation, which at the end modifies the Lane Emden equation. With this simplistic approach, 
we see that the effects with a general equation of state are smaller than those associated to the cosmological constant. In 
particular, the equation of state ui^ — —1/3 displays a null effects since it implies rj — (note that this equation of state 
can also resemble the curvature term in evolution equation for the background). On the other hand, phantom models of dark 
energy, which are associated to equations of state Wx < —1( Caldwell 200^ : Noiiri 20051 ). have quite a strong effect. In Fig. [3] 
we show numerical solutions of Lane- Emden equations for a background dominated with dark energy with u)^ = —2/3 and a 
phantom dark energy with uj^ — —2, with — 10~^. These curves are to be compared with those at fig[ll Clearly equations 
of state with Wx < — 1 will generate larger radius than the case described in the main text. Furthermore, the asymptotic 
behavior of the ratio between the density and pvac is p/pvac It?]- 



3.3 Stability criteria with cosmological constant 



Stability criteria for polytropic configurations can be derived from the virial equation. Using equations (|23|) . (I24|l and (|25p we 
can write the virial theorem (IIOII in terms of the radius R and the mass M: 



2^+2^^ 



)kPc" Mf2 + -TrpvacAfi?^ (5/1 



1) =0, 



(30) 



where we have assumed that the only contribution to the kinetic energy comes from the pressure in the form of /C = |n. 
Note that for pvac = and finite mass, one obtains R cc (5 — n)~^ while for pvac 7^ we would obtain a cubic equation for the 
radius. Instead of solving for the virial radius, we solve for the mass as a function of central density with the help of Eq (|2ip. 
We have 



M = gp," 



«/o/2(5 - n) 



(31) 



The explicit dependence of the mass with respect to the central density splits into two parts: on one hand it has the same form 
as the usual case with A = 0, that is, p'.'^""'/-^"- on the other hand the function Q has a complicated dependence on the central 
density because of the term ^c. With the help of (|2H) and (|3ip we can write a mass-radius relation and the radius-central 
density relation 

M = (^gii^) rI^ , (e^/o)pl^^^ 

Following the stability theorem (see for instance in l| Weinberg 19721 ) 1. the stability criteria can be determined from the 
variations of the mass in equilibrium with respect to the central density. We derive from Eq. (|32p : 

dM 

dpc 



3 / 4\ _i dg 



4(7-1) 



(32) 



(33) 



Stability (instability) stands for dM/dpc > {dM/dpc < 0). This yields a critical value of the polytropic exponent 7crit when 
dM/dpc — given by 



.,(C) = i + 2^ 
3 3 a In pc 



(34) 



in the sense that polytropic configurations are stable under small radial perturbations if 7 > 7crit- It is clear that the second 
term in (|34p also depends on the polytropic index and therefore this equation is essentially a transcendental expression for 

7crit ■ 

It is worth mentionin g that by including the corr ections due to general relativity, the critical value for 7crit is also modified 
as 7crit = (4/3) -I- Rs/R ((Shapiro fc Teukolskv 1983 ) and hence for compact objects the correction to the critical polytropic 
index is stronger from the effects of general relativity than from the effects of the background. Thi s is as we would expect 
it. Stabihty o f relativistic configurations with non-zero cosmological constant has been explored in ijBohmer fc Harko 20051 : 
Bohmer 20oi ). 

Going back to equation (|3ip . we can write the mass of the configuration as M — aAiAf(O), where A/(0) is the mass when 
A = and um ~ aM{Cc,n) is the enhancement factor. Both quantities can be calculated to give 
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n 


Cc = 0.1 


Cc = 0.05 


Cc = 0.001 


1 


(1.12, 1.29) 


(1.05, 1.12) 


(1.001, 1.002) 


1.5 


(',-) 


(1.11, 1.17 


(1.002, 1.003) 


3 


(-,-) 


(-,-) 


(1.022, 1.01) 



Table 1. Numerical values for the enhancement factors (oii,aM) (values for 
llBalaguera-Antolmez fc Nowakowski 2005h ). The symbol — indicates a non defined radius. 



1 have been taken from 



M(0) = (^(5 - n)/("\f("' 



CtM 



Jo J2 



<c/^(5/i - 1)) 



(35) 



where fl 



(n) 



/i(Cc = 0, n) are numerical factors (tabulated in table 1) that can be determined in a straightforward way. 



Similarly, by using Eq (I21|l . the radius can be written as i? = a_Ri?(0), where 



i?(0) = {k{^ 



Jo PC 



Or 



fo 



Jo 



(n) 



/0/2 



/r/r (i-Kc/o^ (5/1-1)) 



(36) 



In table [T] we show the values of the enhancement factors a At and an for different values of (c and different polytropic 
index n. We also show the values of the critical ratio (^crit which separates the configurations with definite ratio such that a 
zero ^1 exist provided that ^c < Ccrit- We will show some examples where the enhancement factors may be relevant in section 
5. 



4 THE ISOTHERMAL SPHERE 

The isothermal sphere is a popular model in astrophysics, either to model large astrophysical and cosmological objects 



1 ne isothermal sphere is a popular model m astropUysics, either to model large astrophysical and cosmological objects 
(galaxies, galaxy clusters) jLvnden-Bell fc Wood 196a : Pensi on 1969^: lYabushita 19681 : Sommer-Larsen.Vedel fc Hellsten 1996 



Chavanis 2001 



Sussmai^_&_JIernan^ez_20M the so-called gravothermal catastrophe (jBinnev fc Tremaine_2987 



Nataraian fc Lvnden-Bell 1997 : Lombardi fc Berti 2001 ) and finally to compare observations with model predictions ( Rines et al. 20o3 : 
Sussmanfc'^emand^ In the limit n ^ oo in the polytropic equation of state one obtains the description for an isother- 

mal sphere, (ideal gas configuration) with 



2 

: a p 



(37) 



where a is the velocity dispersion (cr^ cx T). The pattern we found in section 3 for n ^ 5 gets confirmed here: no finite radius 
of the configurations is found with A, the asymptotic behavior is not p — » as r oo (but rather p — > 2pvac) and, as we will 
show below, other attempts to define a proper finite radius are not satisfactory. 

The results for a finite value of the index n axe defined in the limit n — > oo only asymptotically in the case A = 0. We 
consider this as an acceptable behavior of the density. Because of the limiting case n ^ oo the analysis for the isothermal 
sphere must be done in a slightly different way. As was done in Eq. (|13p . we can integrate the equilibrium equations (Euler 
and Poisson's equations) and obtain an explicit dependence of the density with pvac as 



p(r) = pc exp 



-!^(3>grav(r) - $g.av(0)) 





r 8 2] 


exp 


^vrpvacr 

.OCT 



The resulting differential equation for the density with cosmological constant can be written as 



a^_d_ / 

^2 ^j. y 



) d In p 
dr 



= — 47rp -I- SttPv 



(38) 



(39) 



This differential equation could be treated in the same way as we did for the polytropic equation of state, i.e, by defining 
a new function j ) ~ p/pc, but here we can already use the fact that the cosmological constant introduces scales of density, 
length and time ( Balaguera-Antoh'nez.Bohmer fc Nowakowski 200 6lJ). Let us then define the function V'(0 = lii(p('"oC)/Pvac) 
and r = ro^, with ro the associated length scale. Since we are now scaling the density with pvac, the associated length scale 
ro should be also scaled by the length scale imposed by A: 



ro 



OTA 



13.34 



103 km/s 



Mpc. 



(40) 
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Figure 4. Scaled density p/pvac = e.'^ for the isothermal sphere at different values of central density. The solutions oscillate around the 
value p = 2pvac- 



For an hydrogen cloud with a ~ 4 km/s we have tq ~ 40 kpc which is approximately the radius of an elliptical (EO) galaxy. 
In terms of the function '(/'(Oi ^he differential equation governing the density profile is then written as0 



1 d 



:d^ 



^2 dc y 2 ' 

so that according to Eq. (|38|) we may write 
ip{r/ro) = In ( ) - (cE>grav(r-) - ^Trpv 

V Pvac /V 6 



(41) 



(42) 

From this we can derive different solutions ip depending on the initial condition ^(^ = 0) = ln(pc/pvac) and d^(^)/d^ = at 
^ = 0. In Fig. |4]we show numerical results for the solutions of equation (|41|) using different values values of Pc/pvac- As it is 
the case for n > 5, the radius cannot be defined by searching the first zero of the density i.e. the value such that e''"'-^^^ = 
(including %p — > — oo). In this case, the behavior of the derivative of the density profile changes as compared with the the A = 
case since with increasing ^ the density starts oscillating around the value p = 2pvac such that for ^ ^ oo one has a solution 
p — > 2pvac. This can be checked from (|4ip which corresponds to the first non trivial solution for p. This behavior implies that 
there exist a value of ^ = where the derivative changes sign and hence the validity of the physical condition required for any 
realistic model i.e. dp/dr < should be given up unless we define the size of the configuration as the radius at the value of the 
first local minimum. We will come back to this option below to show that it is not acceptable. A second option would be to 
set the radius at the position where the density acquires for the first time its asymptotic value 2pvac. We could motivate such 
a definition by demanding that the density at the boundary goes smoothly to the background density. This is for tow reasons, 
however, not justified. First, we recall that a positive cosmological constant leads to a repulsive 'force' as it accelerates the 
expansion. A negative cosmological constant could be modeled in a Newtonian sense by a constant positive density which, 
however, is strictly speaking still not a background density. Secondly, if we include the background density p(, we would have 
started with p + pb (with a dynamical equation for p being pb ^constant) in which case the boundary condition would again 
be p{R) = to define the extension of the body (or at least, p ^ oo as r ^ oo). Hence, this second option can be excluded 
on general grounds. In any case as can be seen from Figure|4]both definitions would yield two different values of radius. Since 
the first candidate to define a radius is based upon a physical condition of the configuration, we could expect this definition 
as the more suitable one. However such a definition must be in agreement with the observed values for masses and radius of 



^ Compare with Eq. 374 of l| Chandrasekhar 1967t ) or Eq.l of l|Nataraian , 
density. The factor 1 on the r.h.s of 1411 1 is due to pvac- 



Lvnden-Bell 19971 ) where the density is scaled by the central 
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R/kpc 




Pc/p{R) 


M/Mq 


P (gr/cm^) 


10 


0.00249 


1.0005 


2.17 X 10* 


3.4 X 10-25 


50 


0.01248 


1.0129 


2.7 X 10^0 


1.94 X 10-25 


100 


0.0249 


1.052 


2.11 X 10" 


7.2 X 10^26 



Table 2. Values for {i, pc/p{R) and the mass for different values of radius and for pc = lO^pvac with cr = 300km/s. 

specific configurations and the validity of this definition can be put to test by the total mass of the configuration, given as 

/(Ci) = ree^dC (43) 



103 km/s^ 
Combining (|40|) and (|43|) we can write 



R 



3 



M = 653.55 i^-^j Cr',f(Ci) Mq. (44) 

If we define the radius at the first minimum (see Fig. we find M ^ 2 x 10^ {R/kpcf Mq. Although this might set the 
right order of magnitude for the mass of a EO galaxy if we insist on realistic values for the respective radius, say _R ~ 10 
kpc, the picture changes again as the radius is fixed by (|40p which gives R ~ 5.3 x 10® (cr/lO'^km/s) kpc. In order to get a 
radius of the order of kpc with masses of the order of 10^° we would require a ~ 10-^km/s. This differs by almost eight 
orders of magnitude with the measured va lues for the velocity d ispersion a in elliptical galaxies {a ~ 300 km/s) or with the 



Faber- Jackson Law for velocity dispersion ( Padmananbhan 1993l l. We conclude that defining the radius by the position of the 



first minimum is not a realistic solution. As the last option to get realistic values for the parameters of the configuration we 
consider the brute force method to simply fix the value of ^i. It is understood, however, that this method is not acceptable if 
we insist that the model under consideration has some appealing features (without such features almost any model would be 
phenomenologically viable). Therefore we discuss this option only for completeness. For a configuration with pc = lO^pvac, say 
an elliptical galaxy, we fix the radius at _R ~ 50 kpc in Ea. (|40|l and using a typical value for the velocity dispersion a ~ 300 
km/s we get ~ 0.012 which implies /(^i) ~ 0.036. The mass in Eq. H43|) is then given as M ~ 1.5 X 1O"M0 while the 
density at the boundary is p_R ~ 36000pvac, that is, pc ~ 2.35pR. In table [2] we perform the same exercise for other radii. The 
resulting mean density is in accordance with the observed values of the mean density of astrophysical objects ranging from 
an small elliptical galaxy to a galaxy cluster. However, as mentioned above, the model introduces an arbitrary cut-off and 
cannot be considered as a consistent model of hydrostatic equilibrium. 

In summary, the attempts to define a finite radius for the isothermal sphere fail in the presence of a cosmological constant 
either because such a model fails to reproduce certain phenomenological values (if the definition of the radius is fixed by the 
first minimum) or because the definition is technically speaking quite artificial to the extent of introducing arbitrary cut-offs. 
Note that this conclusion is valid almost for any object as the density of the isothermal sphere with A has a minimum whatever 
the central density we choose. 

5 EXOTIC ASTROPHYSICAL CONFIGURATIONS 

In this section we will probe into the possibilities of exotic, low density configurations. The global interest in such structures 
is twofold. First, the cosmological constant will affect the properties of low density objects. Secondly, A plays effectively the 
role of an external, repulsive force. Hence a relevant issue that arises in this context is to see whether the vacuum energy 
density can partly replace the pressure which essentially is encoded in the parameter n {p — np''). By the word replace we 
mean that we want to explore the possibility of a finite radius as long as the pressure effects are small in the presence of pvac 

5.1 Minimal density configurations 

As mentioned above the effect of a positive cosmological constant on matter is best understood as an external repulsive force. 
In previous sections we have probed into one extreme which describes the situation where a relatively low density object is 
pulled apart by this force (to an extent that we concluded that the isothermal sphere is not a viable model in the presence of 
A). Limiting conditions when this happens were derived. On the other hand, approaching with our parameters these limiting 
conditions, but remaining still on the side of equilibrium, means that relatively low density objects can be still in equilibrium 
thanks to the positive cosmological constant. The best way to investigate low density structures is to to use the lowest possible 
central density. As explained in section 3, for every n there exist a An such that pc J? AnPvmc which defines the lowest central 
density. Certainly, a question of interest is to see what such objects would look like. We start with the parameters of the 



Probing Dark Energy at Astrophysical Scales 13 




configuration. The radius at the critical value ^crit is given by Ea. (|21|l after taking the limit pc AuPvac (see ((28])). It is given 
by 



M 



1/3 



J?crit = 2.175 [^^^i^j /o(CcHt; n)C^ Mpc. (45) 

From fig.[5]we can see that the product /o (Ccrit; M)Ccrit changes a little round the value ~ 0.2 as we change the index n, so that 
these polytropic configurations will have roughly the same radius (for a given mass). This implies that such configurations 
have approximately the same average density p(Ccrit) ~ 1.64/9crit ~ 2.34pvac. That is, such configurations have a mean density 
of the order of the of the critical density of the universe. Given such a density we would, at the first glance, suspect that the 
object described by this density cannot be in equilibrium. However, our result follows strictly from hydrostatic equilibrium 
and therefore there is no doubt that suc h obj ect can theoretically exist. Furthermore, p satisfies the inequalities derived in 
i Nowakowski &: Sanabria fc Garcia 20021 ) and ( Balaguera-Antohnez.Bohmer fc Nowakowski 2005a ) from virial equations and 



Buchdahl inequalities which guarantee that the object is in equilibrium (p > 2pvac). Interestingly, the central density for such 
objects has to be much higher than p as, e.g., for n = 3 we have « 300 and therefore pc > 300pvac. Note that these values 
have been given from the solution of the Lane-Emden equation, which is a consequence of dynamical equations reduced to 
describe our system in a steady state. However, we have not tried to solve explicitly quantities from the virial theorem. This 
makes sense as for Dark Matter Halos (DMH) the parametrized density profiles go often only asymptotically to zero and 
the radius of DMH is defined as a virial radius where the density is approximately two hundred times over the critical one. 
Therefore, an a nalysis using virial equations seems to be adequate here. For constant density Ea. (|30p can be expressed as a 



cubic equation (IBalaguera-Antoh'nez.Bohmer fc Nowakowski 2006bl ) for the radius at which the virial theorem is satisfied (let 



us ignore for these analysis any surface terms coming from the tensor virial equation). However, if the density is not constant, 
this expression becomes a transcendental equation for the dimensionless radius ^vir = Rvii/o.- This equation is 

£2. = 6(5 - n)/g 

(n + l)[3-2<e(5/i-l)/3]' ^^''> 

understanding the functions fi now as integrals up to the value £vir. Once we fix (^crit for a given index n, we use as a first 
guess for the iteration process the value ^i((^crit). In fig[S]we show the behavior of the functions fi and the solutions of Eq. (|46p . 
For these values, Eq. (|45|) gives for n = 3 a radius 

M ^ 



i?vi. = 25.8 ( — ) pc, (47) 



which can be compared with the radius-mass relation derived in the top-hat sphericall collapse ( Padmananbhan 19931 ) 
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where /i is the dimensionless Hubble parameter and z^ir is the redshift of virialization. The resulting average density is then 
of the order of the value predicted by the top-hat sphericall model: 

p = (^rr^) p^^'' ^ 2oopcHt, (49) 

\2TVQcritfo J 

Note with the help of Ea. (|16p and (|45p that the mass can be written as proportional to the parameter k^^^ (introduced in 
the polytropic equation of state Ea. (|ll|l V Therefore k — > is equivalent to choosing a small pressure and, at the same time, 
a small mass which, in case of a relatively small radius, amounts to a diluted configuration with small density and pressure. 
Without A such configurations would be hardly in equilibrium. Hence, for the configuration which has the extension of pc, 
Ea. (|47(l . with one solar mass we conclude that the equilibrium is not fully due to the pressure, but partially maintained also 
by A. This is possible, as A exerts an outwardly directed non local force on the body. Other mean densities, also independent 
of M axe for n = 1.5 and n = 4, respectively: 

p = 15.3pcrit, p = 2.6 X lOVcrit, (50) 

The first value is close to pcrit and therefore also to pvac Certainly, if in this example we choos e a small mass, equivalent to 



choosing a negligible pressure, part of the equilibrium is maintained by the repulsive force of A. In ijBalaguera-Antoh'nez.Bohmer fc Nowakowski 
we found a simple solution of the hydrostatic equation which has a constant density of the order of pvac- The above is a non- 
constant and non-trivial generalization of this solution. 



5.2 Cold white dwarfs 

The neutrino stars which we will discuss in the subsequent subsection are modeled in close analogy to white dwarfs. Therefore 
it makes sense to recall some part of the physics of white dwarfs. In addition we can contrast the example of white dwarfs to 
the low density cases affected by A. 

In the limit where the thermal energy ksT of a (Newtonian) white dwarf is much smaller than the energy at rest of 
the electrons {pp, these configurations can be treated as polytropic configuration with n = 3. This is the ultra-relativistic 
limit where the mass of elec trons is much smaller than Fermi's moment um pp. In the opposite case we obtain a polytrope or 
configuration with n = 3/2. j Weinberg 19721 : Ishapiro fc Teukolskv 19831) . In both cases, the parameter k„ from the polytropic 



equation of state is given as 

where ?n„ is the nucleon mass, me is the electron mass and fj, is the number of nucleons per electron. Using the Newtonian 
limit with cosmological constant, we can derive the mass and radius of these configurations in equilibrium. In the first case, 
for n = 3 the mass is written using pifl as Mo = Q{n = 3) which corresponds approximately to the Chandrasekhar's 
limit (strictly speaking a configuration would have the critical mass, i.e, the Chandrasekhar's limit, if its polytropic index 
7 is such that 7 = 7crit). For this situation one has Mo{n = 3) = 5.87p~^ Mq and Ro{n = 3) = 6.8{pq/ pc)^^^ Rq. 
On the other hand, for n — 3/2 one obtains Mo(n = 3/2) = 3.3 x 10~'^ {pq/ pc)~^^^ p~^^^ Mq and the radius is given by 
Ro{n = 3/2) = Q.27{pQ / pcY^^ p.'e Rq where p© is the mean density of the sun. Since for these configurations the ratio i^c 
is much smaller than 10"'' we see from Fig. [5] that the effects of A are almost negligible. The critical value of the ratio 
gives for n = 3 the inequality pc > 307.69pvac and for n = 3/2 the same limit reads pc > 24.24pvac. Central densities of white 
dwarfs are of the order of lO^gr/cm^ which corresponds to a deviation of nearly thirty orders of magnitude of pvac. 



5.3 Neutrino stars 



An interesting possibility is to determine the eff ects of Pvac on c onfigurat ions formed by light fermions. Su ch config urations 
can be used, for instance, to model galactic halos ( Dolgov fc Han sen 2002; Lattanzi.Ruffini fc Vereshchagin 2 003: Jctz er 19961 : 



Borner 20041 ). While discussing the phenomenological interest of fermion stars below, we intend to describe such a halo. Clearly, 
these kind of systems will maintain equilibrium by counterbalancing gravity with the degeneracy pressure as in a white dwarf. 
For stable configurations, i.e, n — 3/2, one must replace the mass of the electron and nucleon by the mass of the considered 



fermion and set p = 1 in (|51|l and (I3ip . We then get for the mass and the radius: 

Mo = 3.28 X 10^« f^) ' Me = 3.64 x Mq, (52) 
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Figure 6. Mass and radius for difTcrent central densities, for index n - 
eV. 



: 3 and n = 3/2. The masses range from mf = I eV up to my =5 



Ro = 1.31 X 10-^ ( ^ I ( ^ I Mpc = 6.16Cc'" { ' Mpc. 



Pc \ " / eV 
,P0/ \m.f 
On the other hand, for n — 3 one has 



Mo 



5.16 X 10' 



eV 
mf 



Mp 



(53) 



Ro = 0.14 — - 



eV 
ruf 



Pq J 



pc = 



3.1xlO^C^/3(^)^^%c 
Vmf / 



The cases represent, among other, possible cosmological configurations when the fermion mass if of the order of eV, for 
instance, massive neutrinos. In figure |6] some representative values for the choice m/ = 1 eV up to m/ = 5 eV are given. 
Obviously in the case n = 1.5 the values for the mass and radius are sensitive to the choice of m/. Indeed, the dependence on 
the fermion mass is much stronger than on the central density. It is justified to speculate that a relative low density objects, 
affected by A, might exist. If then, as in an example we choose m/ ~ 5 eV and pc ~ 40pvac then the mass comes out as 
10^^ solar masses with a radius of the order of magnitude of half Mpc which might indeed be the dark matter halo of a 
galaxy (or at least part of the halo). As pointed out before, such configurations must have a central density greater than 
24.4pvac in order to be in equilibrium. From table [1] we see that the effect on a configuration with ~ 0.05 is represented 
in an increase in the mass by 17% with respect to A/o and an increment of 11% in the radius. Then the conclusion would be 
that A affects such a dark matter halo. This is to be taken with some caution as the fermions in such a configuration would 
be essentially non-relativistic. Note also that it is not clear if neutrinos make up a large fraction of the halo, however, we 
can also speculate about a low density clustering around luminous matter. Of course, allowing larger central density might 
change the picture. However, the emerging scenario would not necessarily be a viable phenomenological model. For instance, 
changing the value of nif from eV to keV (MeV) would reduce the mass by twelve (twenty four) orders of magnitude which 
is definitely too small to be of interest. We could counterbalance this by increasing the central density by twenty four ( forty 
eight!) orders of magnitude. Such a 'countermeasure' would, however, result in a reduction of the radius by eight (sixteen!) 
orders of magnitude, again a too small length scale to be of importance for dark matter halos. In other words, the example 
with a fermion mass of the order of one eV and low central density is certainly of some phenomenological interest. 

The case n = 3 is similarly stringent. A neutrino mass of 1 — 5 eV gives a mass for the entire object of the order of ten 
to the eighteen solar masses which is too large. A fermion mass of several keV would be suitable for a galaxy halo (a mass of 
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the order Me V and higher would give a too small total mass) . With a relative low central density as before (see Figure |6]) we 
then obtain the right order of magnitude for the halo. But then we will have to live with the fact that such halo reaches up 
to the next large galaxy. Briefly, we touch upon the other possible application of fermions stars which have been discussed as 
candidates for the central object in our galaxy. If we allow the extension of this object to be 120 AU and the mass roughly 
2.6 million solar masses, then the fermion mass would come out as 10^ eV for n — 1.5 (10^ eV for n = 3) and the central 
density as lO^^pvac (lO^^pvac for n = 3). 



5.4 Boson stars 



We end the section by putting forward a speculative question in connection with boson stars. The latter are general relativistic 
geons and can be treated exactly only in general relativistic framework i.e. these kind of configurations are based on the 
interaction of a massive scalar field and gravitation w hich leads to gravitational bounded systems. These objects have been 
also widely discussed as candidates for dark matter i Rufflni fc Bonazzola 19691 : Lai 2004 ). On the other hand, variational 



methods in the connection with the Thomas-Fermi equations give relatively good results even without invoking the whole 
general relativistic formalism. By including A we essentially introduce into the theory a new scale, say in this case a length 
scale ta = l/y/A. The basic parameters of dimension length in a theory with a boson mass m_B and a total mass M — NsmB 
where Nb is the number of bosons are (for a better distinction of the different length scales, we restore in this subsection the 
value of Gjv) 

Li = = GnM, L2=rB = — , L3 = ta = (54) 

niB VA 

The resulting radius of the object's extension can be a combination of these scales i.e. 
R^P oc L„ oc N]^L, 

Ri oc (L?L,)'/^ i?2 oc iVSi?i (55) 
and similar comb ination of higher order. Which one of the combination gets chosen, depends on the details of the model. In a 



close analogy to ( Spruch 1991 ; Eckehard fc Schunck 199^ we can examine this taking into account the presence of a positive 



cosmological constant by considering the energy of such configuration as a two variable function of the mass and the radius 
E = E{R,M) 

E ~ — h -TrGMPvacmBiVB-K (56) 

R H o 

The first term corresponds to the total kinetic energy written as IC — NbP ~ Nb/X and taking A ~ i?. The second term is 
the gravitational potential energy and the third term corresponds to the contribution of the background (see Eq. (H}). By 
treating mass and radius as independent variables (we think this this is the right procedure since the radius will depend on 
the 'external force' due to Pvac), we extremize the energy leading to the following values of mass and radius: 

M ~ ~ 10-- (^) Me, R^( ) ~ 10^ (^) ^ Re (57) 

Such values would lead to a mean density of the o rder of p ^ GP viu: i.e. an extremely low density configuration. The mass given 
in the last expression is the so-called Kaup limit I Spruch 199ll ). Of course, this relative simple treatment does not guarantee 



that the full, general relativistic treatment, will give the same results. Therefore we consider it as a conjecture. However, it is 
also obvious from the the discussion above that low density boson stars are a real possibility worth pursuing with more rigor 
(we intend to do so in the near future). 



5.5 A comparison between ALE profiles and Dark Matter Halos profiles 

It is of some importance to see whether our results from the examination of polytropic hydrostatic equilibrium or from 
the virial equations can be applied to Dark matter configurations. N-body simulations based in a ACDM model of the 
universe show that the density p rofile of virialized Dark Matter Halos (DMH) can be described by a profile of the form 
( Navarro fc Frenk fc White 19961 ) 

p{r) = (2)^-'"p,(r/r,)-™ (1 + r/r,)'"-^ , (58) 

where rs is the characteristic radius (the logarithmic slope is dlnp/dlnrs = —m + ^(m — 3)), ps = p{rs) and the index m 
characterizes the slope of the profile in the central regions of the halo. The mass of the configuration enclosed in the virial radius 
Tvir is given as A/vir ~ 47r(2)'^~'"psrs-F(c) such that one can write ps — (l/3)(2)"'~^AvirpcritC"^-F-"'^(c), where c — r^ii/rs is the 
concentration parameter, Avir is the ratio between the mean density at the time of virialization and the critical density of the 



Probing Dark Energy at Astrophysical Scales 17 




Figure 7. NFW profile (black, solid line) compared to the solutions of ALE equation for different masses and different index n ranging 
from n = 3(blue,dots), n = 4 (red, short-dashed line), n = 4.5 (green, long-dashed line) and n = 4.9 (magenta, dot short-dashed line), in 
the limiting case (Jc = Ccrit • The solution from LE equation is written for the critical value of the parameter fcrit > shown in table [2] The 
vertial black line represents the virial radius from the NFW profile. The vertical colored lines represent the virial radius i?vir (Ccrit) = ^.^vii 
for the different polytropic indices. 



univers e (Ayjr ~ IStt^ in the top-hat model for a flat Einstein-deSitter universe, while Avir ~ 104 in the ACDM cosmological 
model ( Diemand et al 2007 )) and F{c) = x^~"^{l + x)"^~^dx. This universal profile has been widely used in modeling DMH 
in galaxy clusters, and the comparison of these models with a stellar polytropi c-like profile -without the explicit cont ribution 
from the cosmological constant in the Lane-Emden equation- can be found in (jCabral-Roseti et al 2o"oi : lArieli 20oi l. 

Some differences can be described between the ALE and the NFW profiles. First, on fundamental grounds, it is clear that 
the NFW profile does not satisfy the LE equation; a basic reason for this is that dark matter is assumed to be coUisionless 
and is only affected by gravity. However, the fact that the real nature of dark matter is still an unsolved issue leaves an open 
door through which one can introduce interaction between dark matter particles leading to a different equation of state (see 
for instance ( Ren et al. 20061 )). On functional forms, one sees that at the central region the difference is abrupt, since the 
Lane-Emden equation has a flat density profile at r = 0, while the NFW profile has a cuspy profile of the form p ~ r~™. It has 
been widely discussed how such cuspy profile is inconsistent with data showing central regions of clusters with homogeneous 
cores (see discussion at the end of this section). Such small slope in the inner regions can be reproduced by the universal 
profile for the case m = 0, and for the profile derived from the ALE equation (which is just a consequence of initial conditions, 
and hence, its independent of A). In this case, the effects of the cosmological constant can be reduced to explore the outer 
regions of the halos and compare, for instance, the slope of the profiles and the virial radius predicted by each one. 

The density profiles predicted by the ALE equation and the NFW profile (for m — 1) are presented in figlT] for the 
virial radius given by ^vir and in fig |8] for the radius given by ^i, both with ~ Ccrit- Also, the behavior of the radius 
for different values of mass and polytropic indices are given in Fig |9] (with ~ 25.3 (Mvir/IO^'^M©) ' kpc and rvir = 

1A98A~.^/^ (Mvir/lO^^Mo)^^^ Mpc « 255 (Afvir/lO^^Me)^''^ kpc (jCentile fc Tonini fc Salucci 200?! )). 

We see that the virial radius given by the ALE equation is surprisingly close to the virial radius given by the NFW profile 
for n = 3 in the range of masses shown in fig [9] and as pointed in Eq. (|49p . this value yields for this model A^j^^ ~ 200. 
However, as can be seen from the plots, the ALE density profile is almost flat until the virial radius. On the other hand, the 
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m = profile allow us to parameterize it as p(r)/pvac = 2(^(, ^(1 + r/rs) ^, where the concentration parameter c can be written 

In fig llOl and llll we have compared both profiles with the corresponding (ait for polytropic index n = 3 and n = 4.9. For 
the first case we see that the virial radius are of the same order of magnitude than the one predicted by the NFW profile. For 
n = 4.9, these quantities differ by one order of magnitude. In all cases, the slope of the NFW profiles changes faster than the 
ALE profile, which implies that the polytropic configurations enclosed by r^ir display almost a constant density. 

The comparison we made above between the polytropic configuration and the NFW profiles shows that up to the 
cuspy behaviour the polytropic results agree with NFW for the polytropic index n = 3. It is worth pointing out here 
that it is exactly this cuspy behav iour which seems to be at odds with observational facts ( Arieli 20031 : Matos et al. 2005 : 
Hoeft &: Miiecket $ Gottlober 20041 '). Our result for n = 3 is then a good candidate to describe DMH. Indee d, the unde- 
sired feature of t he cuspy behaviour led at least sortie groups to model th e DMH as a polytropic confi guration ("Ari eU 20031: 
Dehnen fc Rose 1993.: Matos et al. 20051: iMcKee 200ll : lOebattista fc Sellwood 1998 : iGonzalez-Casado et al. 2004 : .Yepes et al. 2004) 



Henriksen 20041 : iHoean fc Dalcanton 2000h~ Sometimes it is claimed that a polytropic model is favoured over the results from 
N-body simulation (jZavala et al. 20061 ). However, the oscillatory behavior of the n = 3 solutions of the ALE represents a 
disadvantage when compared with the NFW profiles, although the region of physical interest (below the virial radius) is well 
represented by the solutions of ALE. 



6 CONCLUSIONS 

In this paper we have explored the effects of a positive cosmological constant on the equilibrium and stability of astrophysical 
configurations with a polytropic equation of state. We have found that the radius of these kind of configurations is affected in 
the sense that not all polytropic indices yield configurations with definite radius even in the asymptotic sense. Among other, 
the widely used isothermal sphere model becomes a non- viable model in the presence of A unless we are ready to introduce an 
arbitrary cut-off which renders the model unappealing. Indeed, in this particular case we have tried different definitions of a 
finite radius with the result that none of them seems to be justified, either from the phenomenological or from the theoretical 
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Figure 9. Mass-radius relation for the NFW (red, dashed line) profile compared with the solutions Eq. 1461 (blus, solid line) for four 
different values of n. 



point of view. This is then an interesting global result: A not only afltects quantitatively certain properties of large, low-density 
astrophysical structures, but it also excludes certain commonly used models regardless what density we use. 

For polytropic indexes n < 5 and for certain values of the central density we cannot find a well definite radius. 
These certain values are encoded in a generalization of the equilibrium condition found for spherical configurations (i.e, 
vac) written now as pc > -^npvac- VVe obtain ^\ — 10.8, -^3/2 — 24:.2, — 307.7, ^4 ^ 4000. These values 
set a minimal central density for a given polytropic index n. We have discussed such minimal density configurations and 
determined their average density which strongly depend on n. Interestingly, the radius of such configurations has a connec- 
tion to the length scale which appears in th e Schwarzschild- de Sitter as the maximally possible radius for bound orbits 
( Balaguera-Antolmez.Bohmer fc Nowakowski 2006b. ). In this framework we found also a solution of very low, non-constant 
density. Indeed, the limiting value of the central density in the above equations is a crucial point. Below this value no matter 
can be in equilibrium. However, above this value low density objects can still exist. Both effects are due to A: in the first case 
the external repulsive force is too strong for the matter to be in equilibrium, in the second case this force can counterbalance 
the attractive Newtonian gravity effects (even if the pressure is small). 

Other examples of low density configuration which we examined in some detail are neutrino stars with mass of the order 
1 eV and 1 keV. In there we found that a nowadays dominating cosmological constant affects both the mass as well as the 
radius of such exotic objects. Such effect could change those physical quantities several orders of magnitude. The magnitude 
of such effects, however, depend on the fermionic masses and on the assumption that the fermions in such a configuration 
would be essentially non-relativistic. 

Finally, we made a conjecture regarding boson stars, and used variational methods in connection with the Thomas-Fermi 
equations which could give relatively good results even without invoking the whole general relativistic formalism. We then 
found extremely low density configurations for such astrophysical objects. Notice however, that this conjecture relied purely 
on arguments based on scales and in fact needs a full general relativistic investigation to confirm the results here obtained. 

We have compared polytropic configurations with Dark Matter density profiles from N-body simulations. Surprisingly, 
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we find a reasonable agreement between both approaches for the polytropic index n — 3 and restricting ourselves to the virial 
radius. Our model does not have the the undesired features of the cuspy behaviour of the NFW profiles. 

The importance of the astrophysical properties and configurations found in this article is that they are specific features 
to the existence of a dark energy component. Hence, such configurations (e.g, low density configurations) or properties, if 
ever found in nature would imply a strong evidence for the presence of a dark energy component. Such observations would 
be a completely independent, and so complementary, of other cosmological probes of dark energy such as Supernova la or the 
CMBR. 
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